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We report on a search for direct Kaluza-Klein graviton production in a data sample of 84 pb _1 
of pp collisions at y/s =1.8 TeV, recorded by the Collider Detector at Fermilab. We investigate the 
final state of large missing transverse energy and one or two high energy jets. We compare the data 
with the predictions from a 3 + 1 + n-dimensional Kaluza-Klein scenario in which gravity becomes 
strong at the TeV scale. At 95% confidence level (C.L.) for n = 2, 4, and 6 we exclude an effective 
Planck scale below 1.0, 0.77, and 0.71 TeV, respectively. 

PACS numbers: 14.80.-j, 12.10.-g, 13.85.Rm 



Early attempts to unify gravity and electromag- 
nctism led to the idea of an extra circular spatial 
dimension Because of the periodicity of the ex- 
tra dimension, the metric field of the five-dimensional 
spacetime is Fourier expandable in the extra dimen- 
sion with four-dimensional fields (called Kaluza-Klein 
(KK) modes) as coefficients. 

More recently, Kaluza-Klein theories appear in sce- 
narios of large extra dimensions as introduced by 
Arkani-Hamed, Dimopoulos, and Dvali (ADD) |2j. In 
these theories the standard model gauge theory is con- 
fined to a 3-dimensional domain wall (brane) , embed- 
ded in a higher dimensional compactified bulk space. 
Only gravity propagates in the full bulk space. The 
n compactified extra dimensions are assumed for sim- 
plicity to be "large" circles of common circumference 
R (an n-torus). As a result of compactification, the 
gravitational field that propagates in the bulk can be 
expanded in a series of states known collectively as 
the graviton KK tower. Similar to a particle in a box, 
the momentum of the bulk field is quantized in the 
compactified dimensions. For an observer trapped on 
the brane, each quantum of momentum in the com- 



pactified volume appears as a KK excited state with 
mass m 2 — p n 2 , where p n is the momentum in the 
compactified dimensions, and with identical spin and 
gauge numbers. 

In such a model the Planck scale Mpi , the radius R 
of the compactified space (here assumed to be a torus) , 
and the new effective Planck scale Mry are related 

by i 

M 2 Pl = 8wR n Ml +n , 

where n is the number of extra dimensions. If Mr) 
takes values as low as a few TeV the Higgs naturalness 
problem Q of the standard model can be solved by 
introducing a cutoff not too far above the electroweak 
scale with new physics entering at energies above this 
cutoff. The hierarchy problem of the standard model 
is also recast: the question of why Mpi is so large 
compared to the Z boson mass (Mz) is replaced with 
the question of why R is so large compared to 1/Mz, 
and an ultraviolet hierarchy problem is replaced with 
an infrared one. If we take the most optimistic case of 
M D = 1 TeV and use Mpi ~ 10 19 GeV, we find that 
for n = 1, 2, 4 and 6, R ~ 10 11 m, 1 mm, 10 nm and 
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FIG. 1: Lowest-order Feynman diagrams for the emission 
of real gravitons in pp collisions. 



10 fm, respectively. 

All the states in the KK graviton tower, includ- 
ing the massless state, couple in an identical manner 
with universal strength of Mp^. However there are 
(ER) n massive KK modes that are kinematically ac- 
cessible in a collider process with energy E. The sum 
over the contribution from each KK state removes the 
Planck scale suppression and replaces it by powers of 
the fundamental scale Mo ~ TeV. The interactions 
of the massive KK graviton modes can then be ob- 
served in collider experiments either through their di- 
rect production and emission or through their virtual 
exchange in standard model processes [5j. 

There are three processes in pp collisions that can 
result in the emission of a graviton and a hadronic 
jet: qq — » gG, qg — > qG, and gg — > gG, where q and g 
are quarks and gluons and G is the graviton. The 
lowest-order Feynman diagrams for these processes 
are shown in Figure 1. The calculation of graviton 
emission is based on the effective low-energy theory 
that is valid below the scale Mr>. The correspond- 
ing Feynman rules are cataloged in 0, Q • Since the 
graviton passes through the detector without decaying 
or interacting, the experimental signature is missing 
transverse energy {fr) from the emitted graviton and 
a hadronic jet from the outgoing quark or gluon. 

In this Letter we report the results of a search for 
the direct production of KK graviton modes using the 
rate of events with one or two energetic jets and large 
pr at the Collider Detector at Fermilab (CDF). The 
search is based on 84 ± 4 pb _1 of integrated luminos- 
ity recorded with the CDF detector during the 1994-95 
Tevatron run. 

The CDF detector is described in detail in pj. The 
momenta of charged particles are measured in the 



central tracking chamber (CTC), which is positioned 
inside a 1.4 T superconducting solenoidal magnet. 
Outside the magnet, electromagnetic and hadronic 
calorimeters arranged in a projective tower geome- 
try cover the pseudorapidity region \rj\ < 4.2 Q and 
are used to identify jets. Jets are defined as local- 
ized energy depositions in the calorimeters and are 
reconstructed using an iterative clustering algorithm 
with a fixed cone of radius AR = y/ Arj 2 + A(j> 2 = 0.7 
in 7/ — 4> space 9J. The transverse energy of a jet is 
Et = Es'm6, where E is the scalar sum of energy 
deposited in the calorimeter towers within the cone, 
and 9 is the angle formed by the beam-line and the 
cone axis For this analysis, jets are required to 
have E T > 15 GeV. 

The missing transverse energy is defined as 
the negative vector sum of the transverse energy 
in the electromagnetic and hadronic calorimeters, 
$! T = — s'm8i)ni, where Ei is the energy of the 

j-th tower, hi is a transverse unit vector pointing to 
the center of each tower, and di is the polar angle 
of the tower. The sum extends to \r]\ = 3.6. The 
data sample was selected with an online trigger that 
requires pr = Wt\ > 30 GeV. This is a sample dom- 
inated by instrumental backgrounds and by multijct 
events, where the observed missing energy is largely a 
result of jet mismeasurements and detector resolution. 

The two-stage preselection we use to reject beam 
and detector-related backgrounds, beam halo, and 
cosmic ray events is described in [HJ. Events that 
pass the preselection are then required to have only 
one or two jets with Et > 15 GeV, with at least one 
jet within |^| < 1.1. 

We remove events where the missing energy is due 
to energy flow from a jet to an uninstrumented region 
of the detector by requiring that the second highest 
Et jet does not point in 77 to a detector gap if it is 
within 0.5 radians in <f> of the fr direction. We reduce 
the residual mismeasured multijet backgrounds by re- 
quiring that the minimum 8<p between the pr vector 
and any jet in the event (5<j) m in) is greater than 0.3 ra- 
dians and the z position of the event vertex is within 
60 cm of the nominal interaction point. 

To reduce the physics background contribution from 
electroweak processes with leptons in the final state 
(dominated by W(—> lv)) we require that the two 
highest energy jets are not purely electromagnetic 
(by requiring the electromagnetic fraction f em = 
Epm, l EiTn t < 0.9) and the isolated track multiplicity, 
Ntrk U2 i s zer O- For the final sample, we require 
fir > 80 GeV, E T > 80 GeV for the leading jet and 
Et > 30 GeV for the second jet if there is more than 
one jet in the event. By accepting events with an en- 
ergetic second jet we can reliably normalize the back- 
ground predictions from QCD simulation using the 
jet data, control the systematic uncertainty on the 
signal due to initial/final state radiation (ISR/FSR), 
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TABLE I: The data selection path for the fJr plus one or 
two jets search. 



Selection Requirement 


Events Passing 


Pre-Selection 


300945 


1 < N jet < 2 (cone 0.7, E T > 15 GeV) 




M(l or 2) < 1.1 


157035 


2nd jet gap veto 




Sfimin > 0.3 




\Zvertex\ < 60 CHI 


50938 


fem(l), /em (2) < 0.9 


21012 


NtX=0 


16459 


E T {1) > 80 GeV 




If Njet = 2, E T (2) > 30 GeV 


897 


£t>80 GeV 


284 



and interpret the results with a if-factor (the ratio 
of the cross sections at leading-order (LO) and next- 
to-leading-order (NLO), K = unlo/^lo) included in 
the estimated signal cross section. 

The selection requirements and the number of 
events passing at each stage are summarized in Ta- 
ble 1. 

Background events with missing energy and one or 
two jets arise from standard model sources, predom- 
inantly Z{— ► vv)+jets, W(—> lv)+]ets (£ = r, /x, e), 
and residual QCD production. While Z(— > fp)+jets 
produces real fr+jets, W(—> ^V)+jets mimics our 
signal when the lepton is lost or misidentified. To 
estimate the Z+jets and W+jets background levels 
and their uncertainties in the final sample we normal- 
ize PYTHIA [13j Monte Carlo (MC) predictions using 
the observed Z(— > e + e~)+ jets data sample. QCD 
dijet events mimic our signal when one jet is badly 
mis- measured, resulting in large frp. For the QCD 
predictions we use the HERWIG MC program 0] and 
normalize to the high statistics jet data samples us- 
ing well-balanced dijet events. We estimate additional 
backgrounds from ti, single top and diboson produc- 
tion using MC predictions [ljillj], which we normalize 
using the respective theoretical cross section calcula- 
tions for these processes [li). 

The predicted backgrounds from standard model 
processes are summarized in Table 2. Of the 274 total 
events predicted to pass our selection requirements, 
160 are predicted to come from Z{— * w)+jets, 89 
from the combined W{— * ^)+jets electroweak pro- 
cesses, and 22 from QCD production. Because the 
MC predictions have been normalized to high statis- 
tics data samples, the dominant uncertainty on the 
IF+jets and Z+jets predictions is the 4% luminosity 
uncertainty. The QCD prediction has an additional 
14% uncertainty due to jet energy resolution [ll|. We 
observe 284 events in the data. In Figure 2 the pre- 
dicted standard model ftrr distribution is compared 



TABLE II: The predicted number of events in the final 
sample from standard model sources and the number ob- 
served in the data. 



Background Source 


Predicted Events 


Z{— » vv)+]ets 


160.2 ± 11.5 


r^)+jets 


46.6 ± 5.5 


^)+jets 


23.8 ± 5.0 


e^)+jets 


18.1 ± 4.3 


QCD 


21.7 ± 6.7 


ti, single t, dibosons 


3.9 ± 0.3 


Total predicted 


274.1 ± 15.9 


Observed 


284 
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FIG. 2: Comparison between data (points) and standard 
model predictions (boxes) of the tf/r distribution. There 
are 284 data events, to be compared with 274±16 events 
predicted from standard model sources. The distribution 
is plotted with a variable bin size; the contents for bin 
sizes greater than 10 GeV are normalized accordingly. The 
height of the boxes shows the uncertainty on the standard 
model prediction. 

with the distribution we observe in the data. In Fig- 
ure 3 the same comparison is shown for other kine- 
matic distributions. In both figures the data are con- 
sistent with the expected background. An additional 
contribution from graviton production would result in 
a smooth excess over the background in nearly all the 
kinematic distributions, as shown in Figure 4 for frr. 

We use the PYTHIA MC program to generate 
datasets of graviton emission, using the leading-order 
production cross sections calculated in . The signal 
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FIG. 3: Comparison between data (points) and standard 
model predictions (histogram) of the first and second lead- 
ing jet Et, 5(j>min, and Nj e t distributions. 
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FIG. 4: The predicted fhr distribution from standard 
model processes (histogram) and the one from the ex- 
pected graviton signal (for n = 2, Md = 0.6 TeV, and 
a if-factor of 1.0) added to the standard model (hatched). 
The signal appears as a smooth excess over the standard 
model background. The points are the observed data. 



processes are simulated for n = 2, 4, and 6 extra di- 
mensions, and for a range of values of Mjj>. The signal 
efficiency ranges from 2.9% for two extra dimensions 
to 6.4% for six extra dimensions (due to different rel- 
ative weights of the three production processes of Fig- 
ure 1) and is largely independent of Md. The total 
relative systematic uncertainty on the signal efficiency 
is 25%, mostly due to modelling of ISR/FSR (21%), 
jet energy scale (11%), renormalization scale (8%) and 
parton density functions (2%) [j"fj| . 

Using a Monte Carlo technique to convolute the un- 
certainty on the background estimate with the rel- 
ative systematic uncertainty on the signal efficiency, 
the 95% C.L. upper limit on the number of signal 
events is 62. As shown in Figure 5, for K — 1.0 we 
exclude an effective Planck scale less than 1.00 TeV 
for n = 2, less than 0.77 TeV for n = 4, and less than 
0.71 TeV for n = 6. Recently the DO collaboration 
reported a limit on direct graviton emission using a 
if- factor of 1.3 in the signal cross section [l^. They 
report limits of 0.99 TeV for n = 2, 0.73 TeV for n = 4 
and 0.65 TeV for n — 6. For direct comparison, using 
K = 1.3 our corresponding lower limits on Me> are 
1.06 TeV for n = 2, 0.80 TeV for n = 4, and 0.73 TeV 
for n = 6. 

These are the best limits to date on direct graviton 
emission 0, 0] from the Tevatron. 

Assuming compactification on a torus, the limits 
on Mr, with K = 1.0 correspond to limits on the 
compactification radius of R < 0.48 mm for n = 2, 
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FIG. 5: The three curves show the number of expected 
signal events for n = 2, 4, and 6 extra dimensions as a 
function of the effective Planck scale Md for a 7f-factor of 
1.0. The straight line shows the 95% C.L. upper limit on 
the number of signal events. 



7 



R < 0.014 nm for n = 4, and R < 42 fm for n = 6. 
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